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Abstract 

We calculate the fluxes of 511 KeV photons from the Galactic bulge caused by positrons 
produced in the decays of relic particles with masses less than 100 MeV. In particular, we 
tighten the constraints on sterile neutrinos over a large domain of the mass-mixing angle 
parameter space, where the resulting photon flux would signiflcantly exceed the experimen- 
tal data. At the same time, the observed photon fluxes can be easily caused by decaying 
sterile neutrinos in the mass range 1 MeV< msteriie ^ 50 MeV with the cosmological abun- 
dance typically within 10^^ < i^stcriic ^ 10~^, assuming that i^stcriie comes entirely from the 
conversion of active neutrinos in the early Universe. Other candidates for decaying relics 
such as neutral (pseudo) scalar particles coupled to leptons with the gravitational strength 
can be compatible with the photon flux, and can constitute the main component of cold 
dark matter. 



1 Introduction 



Recent observations of 511 KcV photons from the galactic center produced by positron 
annihilation performed by the SPI spectrometer on the space observatory INTEGRAL [1] 
confirm previously reported photon fiuxes [2] and improve the angular resolution. New data 
are inconsistent with a single point source, although more data with better angular resolu- 
tion are needed to exclude or resolve multiple point sources. An apparent de-localization 
of the source around the galactic center combined with the absence of well-recognized as- 
trophysical sources of positrons that can fill the galactic bulge have prompted speculations 
about a non-Standard Model origin of galactic positrons. 

An interesting explanation was suggested recently in Ref. [3], where galactic positrons 
were attributed to annihilating dark matter. In order to have photon fluxes at the observable 
level, one needs rather large number densities of dark matter particles near the center that 
requires smaller masses than the IGeV - 1 TeV range preferred by most models of weakly 
interacting massive particles. Combined with the additional requirement that the produced 
positrons slow down to non-relativistic velocities before annihilation, this puts dark matter 
particles in a mass range of 1 MeV to 100 MeV [3]. According to the scenario of Ref. [3], the 
galactic positrons can be effectly trapped by the magnetic field and dissipate their energy 
by multiple collisions with the gas. The stopping distance is estimated to be comparable 
or shorter than the "free path" before the annihilation for positrons with energies less 
than about 100 MeV. An interesting spin-off of this scenario is the possibility to probe the 
steepness of the dark matter halo profile near the center of the galaxy. Using the power-like 
parametrization (see e.g. [4]) of dark matter energy density, pdm ~ f^'^i near the center, 
the Oxford-Paris group [3] concluded that 0.4 < 7 < 0.8 is preferred. With more data 
coming in the future with better angular resolution, and with the better understanding of 
conventional (dark- matter unrelated) sources of non-relativistic positrons, the fiux of 511 
KeV photons can provide more detailed information on the distribution of dark matter 
near the center of our galaxy. Even though the explanation of the 511 KeV fine from the 
bulge by annihilating dark matter certainly deserves further considerations, one can easily 
see that this requires very unusual properties of dark matter. Smaller mass and relatively 
large annihilation cross sections to electron-positron pairs would require either extreme fine 
tuning in the models of dark matter (see e.g. [5]) or introduction of new mediating forces 
that would enable more efficient annihilation [6]. 

In this paper, we argue that the decay of quasi-stable relic particles can provide suffi- 
cient amount of positrons to explain the SPI/INTEGRAL data. These relics may be the 
main component of the dark matter, or a sub-dominant fraction of it. From the point of 
view of properties of dark matter, this would not be as restrictive as the scenario with 
annihilating dark matter. At the same time, it is clear that if the decaying dark matter 
can provide sufficient numbers of positrons, the resulting photon angular distribution would 
imply steeper profiles of dark matter near the center than in the scenario with annihilating 
dark matter. 

To substantiate our claims we consider two types of relics that can easily provide a 
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sufficient amount of positrons. Sterile neutrinos in tfie mass range of few KeV - few MeV 
have been extensively studied in a variety of cosmological and astrophysical settings. Their 
cosmological abundance is thoroughly investigated by now as a function of their mass ttt-s, 
mixing angle 9 with active neutrinos [7], as well as their lepton asymmetry and flavor 
[9, 10]. We calculate the decay width of sterile neutrinos and their branching to positrons 
which is always larger than the branching ratio to gamma quanta if the decay into e^e" 
pair is kinematically allowed, >1 MeV. We demonstrate that in the large part of the 
9 — nis parameter space the resulting flux of photons is too large to be compatible with 
observations, and therefore is excluded. Thus, the observed fluxes of 511 KeV photons 
[2] flrmly rule out sterile neutrino dark matter with >1 MeV. On the fringe of this 
parameter space lies an allowed area compatible with the SPI/INTEGRAL observations. 
This area is large enough to span a wide range of the parameter space, compatible with 
sterile neutrino abundances 10~^ < fistcriic ^ 10~^. We also show that the constraints 
on the parameter space of the model from the positron annihilation in the galactic center 
are complimentary to those coming from the excessive cosmological background of gamma 
quanta produced in the decay of sterile neutrinos. 

Other models of unstable relics can be more flexible, and accommodate dark matter 
and observable fluxes of photons at the same time. Here wc consider models with scalar 
and/or pseudoscalar particles with MeV range masses and 0{M^^) couplings to matter. 
Such scalar particles can be motivated by models with extra space-like dimensions, where 
radions may have similar properties. We argue the electron-positron channel can be a 
dominant decay mode for these scalars and that they can be the main component of cold 
dark matter. 

2 Sterile neutrinos as a source of positrons 

Sterile neutrinos are the simplest and the most credible extension of the Standard Model. 
Right-handed neutrinos with a relatively heavy Majorana mass rriR are often used to explain 
the apparent smallness of the mass range for the left-handed neutrinos hinted by various 
neutrino experiments. A priori, the right-handed mass niR can be in a rather large mass 
range, from the GUT scale all the way down to the KeV scale. 

In this paper, we consider a model of a sterile neutrino that is characterized by two 
main parameters, the mass and the mixing angle 9 to the left-handed neutrino: 

N = cos 9Nr + sin 9vl] v = cos 9vl - sin 9Nr, (2.1) 

where refers to a pure SM singlet. We do not confine our discussion to the see-saw model 
and allow 9 and to vary independently, without imposing a my ~ 9'^ma relation. Such 
models were studied extensively in the mass range of few KeV - few MeV [7, 8, 9, 10, 11], 
where sterile neutrinos could constitute entire cold dark matter or at least an appreciable 
portion of it. Note that the left-handed neutrino in (2.1) need not be a mass eigenstate, 
and therefore the sterile neutrino N can have interaction with all three charged leptons. 
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Thus, the decay of N to positrons will have an additional parameter V, \V\ < 1. With this 
parameter, the interaction of sterile neutrinos with a light active neutrino and electrons 
takes the following form: 



_ sin 2^ 



+ higgs terms. 



2 sin 6*14^ cos 6*^1/ V2 V2 

(2.2) 

Note that the existence of interactions with Higgs boson(s) may be important for the 
cosmological production of N only if ms >1 GeV, which is outside of the range of 1 MeV - 
100 MeV that we are interested in here. Thus, we neglect Higgs interactions in any further 
analysis. 

In the mass range of 1 MeV - 100 MeV, the sterile neutrino N can decay into three light 
neutrinos, light neutrino and a e+e" pair, and a neutrino and 7. The decay widths depend 
on the nature of N, which can be both Dirac or Majorana. For a Majorana particle, the 
decay width is twice the decay rate for a Dirac particle. This is simply because two decay 
channels, N — > vuP and N — > i/PP, are allowed for a Majorana neutrino. For the Dirac 
sterile neutrino N the decay width into three light active neutrinos (all three flavors) is 
given by 

^N-*vvi> — To — - — I — , (2-3) 



4 V^V/ 

where Fq is the muon decay width, Fq = (1927r^)~^G|,m^, introduced for convenience. 

The decay of N into one light neutrino and the electron-positron pair may occur due 
to charged and neutral currents. To simplify our formulae, we assume to a few per cent 
accuracy 1 — 4sin^^vK — 0, or sin^ ~ 0.25. In this approximation, and in the limit of 
rrie/ms <S 1, the decay width to ve^e~ summed over all neutrino flavors can be expressed 
as 

sin2 29 
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More accurate formulae that take into account the phase space suppression for ~ 2me 
and more accurate values of sin^ 6^ can be easily derived but not needed for the present 
discussion. The complete result can be found in Ref. [12]. In Eq. (2.4) we take into account 
Z-exchangc and its interference with W -exchange for an arbitrary V . The existence of both 
charged and neutral currents is important, because even for K = 0, T^^uee is different from 
zero, which ensures a non-zero positron width even for "muon" or "tau" sterile neutrinos. 

The decay width into [12, 13] contributes less than 1% to the total decay width, 

r._..r„^(g-^^, 

SO that the total decay width is given by the sum of z/e+e~ and uuu channels, F^r — 
Fjv-n/i/p + Fiv^i^ee- The branching ratio to e^e~ has a very simple form, and depends only 
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on the size of the mixing matrix element V: 



When = 1 the branching ratio to e~^e~ reaches the maximum value of 5/13, while 
V — minimizes it, Br™"- = 1/9. Thus, the intrinsic model dependence of Brg+e- is a 
factor of 3.5. 

Having determined the total decay width of N, we are ready to approximate the cosmo- 
logical abundance of sterile neutrinos by combining the exponential decay factor with the 
power-like pre-exponent reflecting the production of sterile neutrinos in the early Universe. 
Assuming a lepton asymmetry less than 10~^, and using the results of Refs. [7, 10], we 
arrive at 

n.h^ - nos{-r^r} ^ 0.3 (low)'^"^^-^^^^' ^'-^^ 

where r is the age of the Universe. In this formula, JIqs is a fictitious energy density in 
neglection of neutrino decay, introduced here for future convenience. It coincides with an 
actual energy density as long as T^t is small. Note that the pre-exponent in Eq. (2.7) is 
very approximate. In practice, there will be additional dependence on the lepton asymmetry 
of the Universe and on elements of the V matrix [10]. There is an additional hmitation of 
(2.7), namely it works as long as is smaller than T ~ 130 MeV, that is the temperature 
in the early Universe where most of the sterile neutrino production takes place. 

In the mass range between 1 McV and 100 MeV sterile neutrinos behave as cold dark 
matter. Neutrinos will be distributed inside a galaxy with the very same profile as the 
rest of the cold dark matter. Thus, the number density of sterile neutrinos as the function 
of coordinates can be determined as 

ns{r) = ^ , (2.8) 

where Qdm = 0.23 is the "concordance" value for the cold dark matter energy density, and 
PDMi'r) is the cold dark matter energy density in our galaxy. 

Using the above equations, we are ready to estimate the cmissivity of photons that 
originate from the galactic center. Assuming that all positrons quickly slow down in matter 
and annihilate at rest, we estimate the number of emitted 511 KeV photons per volume 
per time, 

A^5ii7(^) - 2n,(r)r^Bre+e-, (2.9) 
from which we can estimate the flux that reaches the Solar system, 

*51l7(^) = J- / N^^^^{T{s))ds, (2.10) 

where ^is the angle from the galactic center and the integration is along the line of sight. 
The total flux is obtained by integrating over the acceptance of the spectrometer. Combin- 
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ing the above equations, one arrives at the following expression: 



. .8.TxlO-3(^)^(^)^4|.r.l) 

xexp{-5x:0-^(^)(^)^4,.r.O)} 

27r J iios VMeV/cmV Vkpc/ 



(2.11) 



Integrals involving the dark matter density are commonly performed with a NFW den- 
sity function [4] which can be simplified to p(r) cxl/r"' for the inner region of the galaxy in 
which we are interested. We have done our calculations for different values of 7. However, 
a better functional form for the NFW density has been developed recently [14] which avoids 
the singularity and thus does a more careful job of dealing with the radial dependence near 
the center. This density function is given by: 



where To — 20/i~^kpc. To fix an overall normalization po we take p — 0.3 GeV/cm^ at 
r = 8.5 kpc, the most commonly used value for the dark matter energy density near the 
Solar system. The values of a which make this expression consistent with previous density 
distributions for larger values of r range between a = 0.1 and a = 0.2. 

The spectrum of the 511 KeV line emission from the galactic center region was obtained 
by the spectrometer SPI on the INTEGRAL gamma-ray observatory. The data suggest an 
azimuthally symmetric galactic bulge component with FWHM of ~ 9° with a 2a uncertainty 
range covering 6°-18°. The flux of 511 KeV photons was measured as $exp — 9.912;! ^ 
ph cm~^ s~^. We have obtained angular distributions and total fluxes expected from the 
decay of sterile neutrinos using the density function (2.12) in the expression for the flux 
(2.11) and averaging over the 2° angular resolution of the spectrometer. By comparing our 
calculated flux to the experimental value we find the effective constraints on sin^ 29 and 
rus. Fig. 1 shows the results of our calculations for ct = 0.1 and 0.2, for the extreme values 
of |K| = 0, 1, with the shaded area providing too large a flux compared to the experimental 
value. Therefore, we can exclude a signiflcant portion of the parameter space in the mass 
interval 1 MeV - 100 MeV regardless what the source of the galactic positrons really is, 
conventional astrophysics or exotic scenarios. We also include in these graphs contours of 
the different values of Qg, given by Eq. (2.7). One can readily see that the flux of 511 KeV 
photons is a very powerful tool for probing sterile neutrinos, as even minuscule abundances 
of sterile neutrinos, fig ~ 10"^, can give observable fluxes. 

To produce a model-independent exclusion plot, we assume the conservative value of 
q; = 0.2 since this value of a yields the lowest count, and we allow |y| to vary. Then the 
conservatively excluded region will be the overlap of the two a — 0.2 plots for the extreme 
values of |y| = 0, 1, shown in Fig. 2. We can also disconnect the result from a particular 




(2.12) 
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Sin2 29 Sin2 29 

Figure 1: Mixing angle - sterile neutrino mass parameter space with the contours of sterile 
neutrino abundance for different halo profiles, a = 0.1 and 0.2, and for different mixing 
angle |V^| = and 1. The dark shaded region covers the combination of parameters that 
creates too large a flux of 511 KeV photons. The boundary of the shaded region gives the 
same flux of photons, as observed by INTEGRAL/SPI. The area bounded by a thick colid 
curve is excluded by the diffuse cosmological background of gamma photons. 
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Figure 2: The exclusion plot for sterile neutrino parameters based on conservative assump- 
tions about the dark matter profile near the center of the Galaxy. The unshaded area 
corresponds to \V\ = 1 and the hatched area to |V^|=0. The solid area is the most conser- 
vative model-independent result, corresponding to a dark matter density which is constant 
inside 3 kpc. 

functional form of the dark-matter density near the center by fixing its value at r = 3 kpc 
and keeping it constant inside that radius while keeping the same parametrization (2.12) 
with a = 0.2 outside of this region. Clearly, this will yield the most conservative outcome, 
and we show it also in Fig. 2. Notice that the assumption about the Majorana nature of 
sterile neutrinos would increase all rates by a factor of 2, and thus lead to a larger excluded 
region. 

Returning to Fig. 1, we observe that the combinations of and 6 that lie at the 
boundary of the shaded region give the photon flux that can explain the INTEGRAL/SPI 
signal. We see that the observed flux can be accommodated with densities in a wide range of 
sterile neutrino abundances, 10~^ < fl < 10~^. Each of these curves contains two branches: 
for a specified value of the mixing angle 9 there are two values of consistent with the 
photon flux. The branch with lower rUs corresponds to the sterile neutrino lifetimes that 
are larger than the age of the Universe, T^r < 1, and low values of fig in this case are 
caused by rather inefficient production of sterile neutrinos in the hot Universe. Such low 
Jig and mixing angles makes this scenario immune to other possible constraints. At the 
same time, for the upper branch of the curve, F^t > 1 and small fig cem also be caused by 
the exponential suppression in (2.7). Consequently, for nis ~ few MeV along this branch 
the energy density of sterile neutrinos at higher redshifts could be of order one compared 
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Figure 3: The angular distribution of 511 KeV photons for two different inner halo profiles, 
a = 0.1 and 0.2, averaged over the 2° resolution of the SPI spectrometer. The experimental 
data were fit with a gaussian of full width at half maximum at a 9° diameter circle (shown 
as a vertical dashed fine), with a 2(7 confidence interval of 6° — 18° (shown as vertical 
dashed-dotted lines). As can be seen, the distribution for a — 0.2 falls well within the 2a 
confidence interval. 



to the rest of matter, and thus is subject to other cosmological constraints such as diffuse 
photon background, (see e.g. [15, 10]). 

If we assume that the fiux of photons seen by INTEGRAL/SPI comes from the decay 
of sterile neutrinos, we can use the angular distribution to determine the preferred cold 
dark matter density profile in the inner galactic region in a similar fashion as in Ref. [3]. 
Averaging over the angular resolution of the spectrometer and using the common flux 
normalization at the zero angle result in the angular profiles shown in Fig. 3 for a = 0.1 
and a — 0.2. The distribution for a — 0.2 falls well within the 2a confidence interval, and 
thus is clearly favored by the data over smaller values of a. The angular distribution that 
we calculated with the power-like density function p{r) ccl/r"' for 7 = 1 comes very close 
to the one for a = 0.2. 

The photon width Tj^i^^^, although smaller than F^y^j/ee, should lead to the additional 
photon line at = ms/2. The intensity of m,/2 line, can nevertheless come close to that 
of 511 KeV if 777.5 is large because of m<i/(2r;ie) enhancement factor. Unfortunately, the 
prediction of photon fiuxes at ms/2 relative to 511 KeV, 

^o.5m,7 _ ^N-^i^-r _ 27q; _ 0.031 

$5117 ~ 2r^_,ee- ~ 7r(4|y|2 + 1) ~ A\V\'' + 1 ' ^ ■ ^ 
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varies by a factor of 5 depending on \V\. 

It is important to compare the limits from positron production in the decays of the 
sterile neutrinos with the hmits derived from a comparison to the cosmological diffuse 
gamma background. Ref. [9] estimates the diffuse photon background produced by the 
radiative sterile neutrinos with an additional assumpiton of sterile neutrinos constituting 
the dominant part of dark matter. As a a first approach, we can effectively compare their 
results by replacing their assumed neutrino abundance of 0.3 by our calculated 10~^ < 
^sterile ^ 10~^, which rclaxcs their constraints on sin^^^ by many orders of magnitude. For 
the diffuse photon background this meant, effectively, a bound on the hfetime r > 4 x 10^^ 
years, but as can be seen from our Eq. 3.19 our sensitivity to r is at the level of 5 x 10^'' 
years asuming the same abundance, 1 MeV neutrinos and a minimum branching ratio 
to e+e" of 0.1. Therefore, it appears that 511 KeV photon fiux provides more stringent 
bounds than the diffuse photon background by one to two orders of magnitude. It is 
nevertheless important to consider this issue in more detail, as the experimental data on 
photon background used in Ref. [9] are clearly outdated. 

The spectral density of the photon flux can be expressed in terms of the sterile neu- 
trino mass rris, the branching ratio to photons, the neutrino lifetime, the "initial" neutrino 
abundance Qqs, the age of the Universe r and its energy density today po, 

= — fiosPoBr^- — VMTeK.^{-Vr,T{2E mf^}. 2.14 

ah Att m \ m J 

The maximum in this spectrum corresponds to photon energies 

_/ m/2 for 3r„r<l, 
^"^'^^ ~ \ (3r„r)-im/2 for 3r„r > 1 ^ ^' 

On the experimental side, [16], we assume a conservative bound on the spectrum of the 
cosmological diffuse gamma ray background, 

^ MeV str-^cm-^s-^ (2.16) 



Plugging in the expression for E^ax (2-15) to the photon spectrum (2.14), we require 
{d^^/dE)max to be smaller than {d^j/dE)exp- This results in the following constraint 
on the parameters of the model, 

{ 9+lfv|2 ^05^^3;exp{-|x} for x < 1, 
(2.17) 
9^^0os/i^a;-3/2exp{-|a;-V2} for x > 1, 

where x — SF^^t is introduced for concision. This constraint simplifies when x <^ 1, in which 
case we obtain a very powerful limit on the combination of masses and mixing angles. 



/ sin^ 26\^ / rris \ 



V 10-19 J \1 MeV 



j ^0.8 (2.18) 
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where V is taken to be one. 

Plotting this constraint in Figure 1, we superimpose it on the hmits obtained from the 
511 KeV hne. We observe that the 511 KeV hne typically provides better sensitivity to 
parameters of the model when T^r < 1, while the cosmological background of gamma 
can provide more restrictive bounds if Tnt > 1. Unlike the limit from 511 KeV photons, 
the cosmological background constraint is not sensitive to the details of the dark matter 
distribution in the galactic center. 



3 Positron flux from decaying dark matter 

In this section, we generalize the analysis of the previous section to an arbitrary decaying 
relic, without immediately specifying its cosmological origin and keeping its abundance fts 
as a free parameter, Qs ^ ^dm- The equality would mean that this relic is the main 
component of the cold dark matter. The flux of the galactic 511 KcV photons, normalized 
to the measured value, can be written as a function of the mass ms, the decay width into 
e'^e" and Q5: 

tK.r...-,a) ^ (5 , ^ r....) ik. (3.19) 



4„p ' " ' ms 0.23 

In this equation, a = 0.2 was used. We observe that even if dark matter particles have much 
longer lifetimes than the age of the Universe, they are still capable of providing significant 
numbers of positrons. 

Sterile neutrinos in the mass range of 1 MeV - 50 MeV can account for the galactic 
positrons but cannot constitute any appreciable fraction of the dark matter if they are 
produced thermally from active neutrino species under the assumption of a small lepton 
asymmetry. To this end, it is interesting to relax the assumptions that led to Eq. (2.7) and 
using Eqs. (3.19) and (2.4) estimate the resulting photon fiux. 



--^-^ + j llOMevJ 0.23- 



- exp 



Combining the analysis of Ref. [10] at \V\ = 1 and the lepton asymmetry L = 0.1 with 
Eq. (3.20), we find that all points on the 9 — rUs parameter space that correspond to 
> 10"^ are firmly ruled out by the excessive photon fiux. Therefore, even for an arbi- 
trary lepton asymmetry, sterile neutrinos with rUs > 1 MeV cannot constitute a significant 
fraction of dark matter. The combination of both features, the sterile neutrino dark mat- 
ter [ttIs > IMeV) and acceptable positron numbers from its decay, may only come at the 
price of abandoning the conventional mechanism of thermal transfer of active neutrinos into 
sterile, and assuming some initial sterile neutrino energy density due to either non-thermal 
production after infiation and/or decay of some heavy non-SM particles into N before the 
electroweak epoch. 
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Another example of an unstable relic that could produce positrons is an extremely 
weakly interacting scalar or pseudoscalar particle (f) in the same MeV mass range. Quite 
generically, we can write its interaction with operators composed from the Standard Model 
fields as 

= ld,<pdy - \m\<^^ + ^^Of ^ (3-21) 

i 

where the leading order dimension 4 operators that induce the decay of are included. In 
models with extra space-like dimensions, a scalar particle called radion would couple to the 

trace of the stress-energy tensor of visible matter. If represents a radion, a scalar mode 
corresponding to the fluctuation of the volume of extra dimension, the coupling constant 
can vary within a large range, from Mpi or even larger scale down to Myy with the latter 
being the case in the Randall-Sundrum model [17] . 

For 1 MeV< ^ 100 MeV, the leading interaction with fermions is given by (f)M~^meee 
for scalars and 0M~^ me 6^756 for pseudoscalars. In both cases, the width to e+e~ is given 

by 

i ^^e+e- - g^^2 - X yrj ^^^y J , [6.22) 

where the threshold effects are neglected, 2mg < m^. Whether Eq. (3.22) provides the 
dominant contribution to the total decay width depends on whether or not the anomalous 
interaction with photons, 4)Fnj^F'^'^ or 0F^,^F''^, is generated at one loop. This depends on 
the UV completion of the effective Lagrangian (3.21), and both outcomes, zero or non-zero 
couplings to photons, are possible (for the scalar case see e.g. [18] and references therein). 
In particular, the one-loop anomalous coupling of scalars and photons is not generated in the 
models where the interaction term of the scalar with matter is obtained from a conformal 
rescahng of Brans-Dicke theory with the mass term for 0. If (pF"^^ arises at two loops or 
not at all, then the total decay width is dominated by V(j,^e+e- ■ For the pseudoscalar case, 
one could also avoid couplings to photons by postulating a derivative interaction with the 
SM flelds, i.e. (9^</)-?/'7^75-0. Taking the point of view that e^e~ represents the main decay 
channel, we combine Eqs. (3.19) and (3.22) to provide an estimate of the photon flux from 
the Galactic bulge, 

m..%)^,, /lO^^GeVy 

$exp \ ) 0.23' ^ ' 

which turns out to be independent of m^. It is remarkable that as large as the Planck 
scale, 

M, ~ 5 X 10^^ GeV^n0/O.23, (3.24) 

can be probed this way for maximal Q^! 

Is it realistic to have ~ 0.23? There are many possible scenarios of non-thermal 
production of scalars and pseudoscalars. In fact, the scalar field energy density can be 
over-produced easily leading to the closure of the Universe (moduli problem). This is, of 
course, very much model-dependent, and can be easily avoided in some scenarios. Thermal 
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production can also be quite efficient if in addition to very weak decay terms the scalar ^ has 
appreciable bilinear couplings to the Standard Model fields. For example, the Lagrangian 
[19, 5, 20], 

= \d^W4> - \{rnl(t>^ - Xcj>'H^H) + ^m,ee, (3.25) 

where H is the Standard Model Higgs doublet, has the potential of providing a dark matter 
candidate 0. To ensure the right abundance of (j) one can choose A to be large, in which 
case is thermalized above and depletes its abundance via annihilation at the freeze- 
out. In this case however, there is a direct experimental bound on that requires it to 
be heavier than 400 MeV [21], thus making it unsuitable for explaining the positron flux. 
Another possibility for obtaining a cosmologically required value of is to assume a small 
value of A, initial absence of scalars, and their thermal production at T ~ My^ [20]. This 
option is indeed viable, and A ~ 10"'^ would lead to ~ 0.2. Thus, we see that the 
scalar field model (3.25) can account for the existence of cold dark matter and lead to the 
observable fluxes if condition (3.24) is satisfied. In the original formulation of this model 
without the decay term the stability of the scalar was ascribed to either Z2 symmetry [5] 
or global U{1) symmetry if is a complex scalar [19]. The appearance of the decay terms 
~ MpiVmeCe might be speculated to arise due to the violation of these symmetries in string 
theory /quantum gravity at the Planck scale. 

4 Discussion and Conclusions 

We have shown that the decaying relics that may or may not constitute the main com- 
ponent of dark matter can easily account for the flux of galactic photons observed by the 
INTEGRAL/SPI instrument [1]. The lifetimes of these relics of order 10^''' yr, i.e. much 
larger than the age of the Universe, can be probed this way. This is also longer than the 
lifetime of unstable ~MeV-scale relics that can be probed with diffuse cosmic gamma-ray 
background [9]. This sensitivity motivates further theoretical studies of the flux of 511 KeV 
photons from the bulge. 

We have studied in detail the model of sterile neutrinos produced in the early Universe 
by conversion from active neutrino species. We have shown that if the decay to e^e~ pairs 
is kinematically allowed, the branching ratio in this mode is never smaller than 10% even 
if the sterile neutrinos have no charged currents with electrons (K = 0). Using previous 
results of abundance calculation [7, 10] and the total decay width of the sterile neutrino 
and its branching to positrons obtained in this work, we were able to calculate the flux 
of the 511 KeV photons as a function of the mixing angle 9 with active species and the 
sterile neutrino mass m^. Over a large domain of the parameter space for 1 MeV< rus ^50 
MeV this flux turns out to be significantly larger than the observed value. This enables 
us to exclude a significant portion of the parameter space of sterile neutrinos regardless of 
whether the source of the galactic positrons is astrophysical or it is due to exotic physics. 
In particular, we are able to rule out any cosmologically significant quantities of sterile 
neutrinos, Q,s < 10~^ with masses larger than 1 MeV. 
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The boundary of the excluded region corresponds to fluxes that are comparable to the 
observed value. If indeed the experimental results can be interpreted as the consequence of 
positrons from the decaying sterile neutrinos, the photon flux points to a sterile neutrino 
abundance in the range of 10~^ <^ ^^stenie ^ 10~^. It is also interesting to note that the 
upper value for the mass of sterile neutrinos (~50 MeV) compatible with the observed flux 
of photons comes naturally into this analysis, and not as an external requirement as in the 
case of annihilating dark matter [3]. The radial distribution of the decaying relics inside 
the galactic halo should trace the distribution of the dark matter. Thus, if indeed future 
data support the link of galactic positrons to exotic physics {e.g. by firmly establishing 
their diffuse origin as opposed to point-like sources), the two scenarios, decaying relics and 
annihilating dark matter, can be distinguished by the angular distribution of the signal. 
The signal from decaying relics is proportional to the number density of the dark matter 
particles n(r) while the annihilating dark matter has n^(r) dependence. Thus, the decaying 
relics would suggest steeper Galaxy profiles than the scenario with the annihilation of dark 
matter. In both scenarios, decaying relics and annihilating dark matter, one would generally 
expect the additional sharp line due to the annihilation or decay to photons. In the case of 
decaying sterile neutrinos, the flux of photons with energies of 0.5ms is predicted to be at 
the level of few percent of the flux of 511 KeV neutrinos. Finally, we would like to comment 
that the scenario with annihilating dark matter can be searched for experimentally using 
the existing data on e+e" colliders. Since the results of Ref. [3] suggest an appreciable low- 
energy dark matter annihilation cross section into e'^e~ , one would naturally look for the 
missing energy signal (plus a hard photon) in the e~^e~ collisions at higher energies where 
the SM model cross sections become smaller and while the cross section of the MeV-mass 
dark matter production would become larger. 
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